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This paper studies the nature of the effective velocity boundary conditions for liquid flow over a 
plane boundary on which small free-slip islands are randomly distributed. It is found that, to lowest 
order in the area fraction j3 covered by free-slip regions with characteristic size a, a macroscopic 
Navier-type slip condition emerges with a slip length of the order of a/3. The study is motivated 
by recent experiments which suggest that gas nano-bubbles may form on solid walls and may be 
responsible for the appearance of a partial slip boundary conditions for liquid flow. The results are 
also relevant for ultra-hydrophobic surfaces exploiting the so-called "lotus effect" . 
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I. INTRODUCTION 



The recent blossoming of research in micro-fluidics has prompted a renewed interest in the 
possibility of slip boundary conditions at the contact of a liquid with a solid wall 0, Q] . While 
many experiments have provided evidence for a violation of the classical no-slip boundary condition 
at small spatial scales y, 0, 0, |28|, |32l H3, OS |3|| the physical mechanisms responsible for this 
phenomenon are still unclear. An interesting possibility is the recent discovery of what appear to 
be small gas nano-bubbles or pockets attached to the wall [3, lJ LLlJ LLil ISJ- f27L kSlj . The evidence 
for the existence of these nano-bubbles is somewhat indirect, but nevertheless compelling. It is also 
hypothesized and, sometimes, experimentally verified (see Watanabe et al. that gas pockets 

may form in cracks or other imperfections of the solid wall, thereby decreasing the overall wall 
stress. 

In order to explore the macroscopic consequences of the existence of such drag-reducing gaseous 
structures on a solid wall, in this study we consider by statistical means the effective velocity 
boundary condition produced by a random distribution of small free-slip regions on an otherwise 
no-slip boundary. We consider both the three-dimensional problem, in which the regions are 
equal disks, and the two-dimensional problem, in which they are strips oriented perpendicularly 
to the flow. While idealized, these geometries provide some insight into the macroscopic effects of 
randomly distributed microscopic free-slip regions. 

We find that, away from the wall, the velocity field appears to satisfy a partial-slip condition 
with a slip length proportional, to leading order, to the product of the length scale a of the free- 
slip islands and the area fraction (3 covered by them. After deriving a general result, we solve the 
problem to first order accuracy in (3 for both the two- and three- dimensional situations. 

As discussed in section IVIII our results are consistent with those of a recent paper by Lauga & 
Stone 01 > wn0 assumed a periodic distribution of free-slip patches on a boundary, as well as those 
of an older paper by Philip wn0 similarly investigated the effect of free-slip strips arranged 
periodically on a plane wall parallel or orthogonal to the direction of the flow. 

The present results are also related to so-called "lotus effect" pj exploited to obtain ultra- 
hydrophobic surfaces. Such surfaces are manufactured by covering a solid boundary with an array 
of hydrophobic micron-size posts which, due to the effect of surface tension, prevent a complete 
wetting of the wall 0, 0, 0] . In the space between the posts the liquid remains suspended away 
from the wall with its surface in contact only with the ambient gas and a concomitant reduction 
in the mean traction per unit area. Another instance of drag reduction by a similar mechanism 
has also been reported in Watanabe et al. [3j| . These authors studied the pressure drop in the 
flow of a viscous liquid in a tube the wall of which contained many fine grooves which prevented 
a complete wetting of the boundary. 

The approach used in this paper is mainly suggested by the theory of multiple scattering 0, 
l29l l3fl ] and was used before to derive the effective boundary conditions at a rough surface for the 
Laplace and Stokes problems [22H23I ] 
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FIG. 1: Solid no-slip boundary with a random distribution of equal circular free-slip areas. 



II. FORMULATION 

We consider the flow in the neighborhood of a locally plane boundary [33] B with a compos- 
ite micro-structure which dictates free-slip conditions on certain areas s 1 , s 2 , . . . , s N and no-slip 
conditions on the remainder B — U^ =1 s Q (figure QJ. If each "island" s a is sufficiently small, and 
U^ =1 s a is also sufficiently small (both in a sense to be made precise later), near the boundary the 
flow is described by the Stokes equations: 

Vp = yN 2 u, V u = 0, (1) 

in which p and u are the pressure and velocity fields and /j, the viscosity. On the free-slip regions 
u satisfies the condition of vanishing tangential stress: 

tj-(T-n) = xes a a = 1, 2, N J = 2, 3 (2) 

where ti and £3 are two unit vectors in the plane and t the viscous stress tensor, while, on the 
rest of the surface, 

u = x^(J% =1 s a . (3) 

The normal velocity vanishes everywhere on B. 
We start by decomposing the solution (p, u) as 

N N 
a—1 a—1 

Here u° and p° are the (deterministic) solution satisfying the usual no slip condition on the entire 
boundary B while the fields (q a ,v a ) account for the effect of the a-th island. We define these 
local fields so that v a vanishes everywhere on B except on s a , where it is such that the free-slip 
condition (J2J is satisfied. To express this condition it is convenient to define 

w a = u° + vP . r a = p° + J2 <f > (5) 

(3=ia 

so that, for every a = 1, 2, . . ., N, 

u = v a + w a , p = q a + r a . (6) 

On s a , then, v a satisfies 

tj ■ {r va ■h) =-tj- (r wa ■ ft) xes a J = 2, 3 (7) 
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where 

T va = n [Vv a + (Vv a f] , r wa = n [Vw a + {Vw a f] , (8) 

the superscript T denoting the transpose. Clearly 

v a ^0, q a ^0 as |a;-y Q |^oo (9) 

with y a a reference point on the a-th island. It is evident that both fields v a and w a satisfy 
the Stokes equations. In the terminology of multiple scattering, they are often referred to as the 
'scattered' and 'incident' fields, respectively pll2ll|. 

III. AVERAGING 

We assume that the free-slip islands are identical circular disks with radius a, centered at y a , with 
a = 1, 2, ... , N. We make use of the method of ensemble averaging and consider an ensemble 
of surfaces diffcrring from each other only in the arrangement of the N free-slip islands. Each 
arrangement, or configuration, is denoted by C N = (y 1 , y 2 , y N ). A particular configuration will 
then occur with a probability P(C N ) = P(N) normalized according to: 

± J dV ■ • ■ / d 2 y N P(y\ ... tV ») = ±.J dC N P(N) = 1 . (10) 

The ensemble-average velocity is defined as 

(«>(*) = ^J dC N P(N)u(x\N) (11) 

where the notation u(x\N) stresses the dependence of the exact field not only on the point x, 
but also on the configuration of the N islands. In view of the fact that u° is independent of the 
configuration of the disks, substitution of the decomposition (@J into (jl l|t gives 

1 N f 

(u)(x) =u°(x) + — J2 dC N P(N)v«(x\N). (12) 

Since the disks are identical, each one gives the same contribution to the integral. Upon introducing 
the conditional probability P(N - l|y x ) defined so that P(N) = P(y 1 ) P(N - 1| y 1 ), we may 
therefore write 

(«>(*) = u "(x) + ^-^JdC N P(y 1 )P(N-l\y 1 )v 1 (x\y\N-l) (13) 
or, in terms of the conditional average 

(v'hWy 1 ) = ^p^yy J dC-ipiN-lly^Wv^N-l), (14) 



(u)(aO = u°(x)+ f d 2 yP{y){v) 1 {x\y) (15) 

JB 

where the integral is over the entire boundary. For convenience, here and in the following, we drop 
the superscript 1 on the quantities referring to disc 1 . Since v a and q a satisfy the Stokes equations 
everywhere, so do (v)i and (q)i- The boundary conditions are 

=0 x i s (16) 

while 



tj- ((t v )! -n) = -tj- {{t w ) x -h) xes J = 2,3 



(17) 
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Note that 

<t$)i = + (^)i) = n(fliWi+4Wi) , (18) 

and similarly for (rj fc )i since averaging and differentiation commute as is evident from the definition 
(|lip. The normal velocity vanishes everywhere: 

= n- =0 a;eB. (19) 

It may be noted that P(x) is just the number density of free-slip islands per unit surface area 
of the boundary; the area fraction (3 covered by these islands is 

P(x) = f P(y)d 2 y ~ na 2 P(x) + 0{a 2 /L 2 ) (20) 

•J \x — y\ <a 

where L, assumed much greater than a, is the characteristic length scale for variations of the 
number density. 

The framework just described can be readily extended to disks of unequal radius, and to non- 
isotropic islands such as ellipses. In both cases the probability density would depend on a suitably 
enlarged list of variables such as the disk radius, the characteristic size, orientation and aspect 
ratio of the ellipses, and so on. 



IV. THE EFFECTIVE BOUNDARY CONDITION 



Now we derive a formal expression for the effective boundary condition on B. To this end, let 
Gjj(y;x) be the Green's tensor for the Stokes problem vanishing at infinity and on the plane 
boundary B. Then 



te>i(x|y) = / [(-qfn + (T° ■n) z ) 1 (s\y)G^(s;x) + (v 3 ) 1 (z\y)T^ k (s- 1 x)n k ]d 2 

JB 



(21) 



where T^ k is the stress Green's function associated to G^ and the integral is extended over the 
entire plane boundary |l2t Il9l | This formula can be considerably simplified recalling that, on the 
boundary, v vanishes everywhere outside s while G w vanishes everywhere. Furthermore, on s, the 
tangential tractions also vanish. Hence, upon taking the xi-axis along the normal with x\ = on 
the plane, we have 



(v^xly) = f(v j ) 1 (s\y)T%[(s;x)d 2 s 

J S 



(22) 



where now the integration is extended only over the free-slip island. We now consider points x 
such that \x — s\ ^> a, but such that \x — s\ is sufficiently small to be in the Stokes region adjacent 
to the boundary. It can be verified that, in this range, we have 



T™(s;x) - 2T ijl (y;x) 



where Tijk is the free-space stress Green's function: 



1 + 



x - s 



Tijk(y,x) 



3 (Vi - x i ){y J - Xj)(y k - x k ) 



-in 



\y - x\ 



Thus, ||22J) becomes 



where 



(vj)i{x\y) ~ 2na 2 T ljl {y-x)V l {y) 



Vi(y) = \ \ (v t ) 1 (s\y)d 2 s 

7ra J\s-y\<a 



(23) 

(24) 

(25) 
(26) 
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is the average velocity over the disk centered at y. Note that V\ = as v\ = 0. This result may 
now be inserted into the expression (|15fl for the average field to find 

(u 3 )(x) = u](x) + 2Tia 2 J #vP(y)T iiX (y;x)Vi(y). (27) 

We now take the 'inner limit' of 127fl by letting the field point x approach B to find (see e.g. 
Pozrikidis [13 pp. 23 and 27) 



1 „ 
2 



lira Tiji (y; x) = -5 l3 5{x - y) (28) 



so that 

(uy)(x) = Tra 2 P(x) V(x) (29) 

where tin is the velocity component parallel to the boundary. Since the problem is linear, a 
dimensionless tensor Wij must exist such that 

Vi = jW ij ((r w ) 1 -n) j (30) 

so that the average field satisfies the partial slip condition 

(tillXsc) = —P(x) W ■ ((Ox ■ n) ~ ^ W ■ «t")i • h) . (31) 

This equation shows that the slip length is of the order of a(3. 

We can now be more specific about the assumption made at the beginning of section |n] as to 
the validity of the Stokes equations near the wall. The condition for this assumption is evidently 
that the Reynolds number 

Re = « (32 ) 

with v the kinematic viscosity, be sufficiently small. Equation l|3()(l shows that |V| is of the order 
of a//i times the magnitude of the wall shear stress; a precise result in a particular case is derived 
in Appendix A. 



V. FIRST-ORDER PROBLEM 



While exact, the result (|31|l expresses the effective boundary condition on the unconditionally 
averaged field (it) in terms of the conditionally averaged wall stress {t w )\. In order to obtain 
the conditionally averaged velocity (tt)i necessary to evaluate this quantity, one would need an 
effective boundary condition which would involve the wall stress averaged conditionally with the 
position of two free-slip islands prescribed, and so on. This is the well-known closure problem that 
arises in ensemble averaging. An explicit solution can only be found by truncating somehow the 
resulting hierarchy of equations. 

The lowest-order non-trivial truncation can be effected with an accuracy of first order in the area 
fraction /3. It is well known that, in this limit, the average 'incident' (w) may be approximated by 
the unconditional average (it) so that 

<«||>(as) = ^\N-[(r)(x)-n}+o(l3). (33) 
A' 

If the density of the islands is small, since w accounts for the effect of all the other islands on the 
one centered at y, (w)i is slowly varying near y so that 

(ti;>i(aO = (w)x(y) + [(x - y) ■ V](w)x{y) + ... (34) 

and, therefore, 



( T jk)i = V(dj(wk)i+dk(wj)i) - fJ>(dj(u k ) +d k (uj)) = {r jk ) 



(35) 



G 



is approximately constant over the island \x — y\ < a. The velocity field (v)i is therefore the 
solution of the Stokes equations vanishing at infinity and whose normal component vanishes 
on the entire plane; the two tangential components vanish for \x — y\ > a while, for J — 2, 3 and 
\x-y\ <a 

tj ■ ([V(«)i + (V(u)i) T ] ■ n) = --tj ■ ((r) ■ h) = const. (36) 

This problem is solved in the Appendix A where it is shown that 

Wij = ^ (37) 

so that the effective boundary condition (|31|l becomes 

(«„)(*) = ^[3(x)((t)(x) ■ h) + o(P) . (38) 

It may be expected that, if the islands had an intrinsic direction (e.g., an elliptical shape) and 
were not randomly oriented, the tensor would not be isotropic so that the average surface 
traction and surface velocity would not be collinear. 



VI. THE TWO-DIMENSIONAL CASE 



The previous analysis can also be applied to the analogous two-dimensional case, i.e. a surface 
with a random distribution of parallel, or nearly parallel, free-shear strips of width a oriented 
perpendicular to the flow direction. The developments at the beginning of section lTVl are still valid 
and we may start from i|22|) noting that, in place of l|24|) . we have 

T ijk (y ; x) = i ^fafe^ • ( 39 ) 
7r \y x\ 

so that (|22|l becomes, in this case, 

(vj)i(x\Z) ~ 2aT i21 (fcx)V 2 (0. (40) 
where ^ is the coordinate in the direction parallel to the plane. Here 



V&) = - I («2>i(Cl0dC (41) 

lc-ei<o 



a 



is again the average velocity over the strip centered at £. The expression l|15|l for the average field 
is modified to 



(«,•) (x) = u oj (x) + 2a J dH P(0 T j21 x)V 2 (0 . (42) 

The analog of l|28|l is still valid so that 

(«2>(0 = aP(®V 2 (0 (43) 

where u 2 is the velocity component parallel to the boundary. 

As before, from the linearity of the problem we deduce the existence of a dimensionless quantity 
W such that 

V = >(r;)i (44) 

so that the average field satisfies the partial slip condition 

(«)(0 = -«0^)i (45) 
H 
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where we have used the fact that the fraction of the boundary covered by the free-slip strips is now 
given by 



m= P(Qd( ~ aP(0+O(a'/L'). (46) 

J\C-£\<a 

The solution of the problem in the dilute limit is given in Appendix B. One find 

W = £ (47) 
so that the effective boundary condition becomes 

(«>(0 =^^P(0(t* v )+o((3) (48) 

VII. CONCLUSIONS 

We have derived an effective velocity boundary condition on a wall covered by a random ar- 
rangement of free-slip disks or two-dimensional strips. For the case of disks we have found that, to 
leading order in the fraction (3 of the unit area covered by the disks, the velocity satisfies a Navier 
partial slip condition with a slip length £ given by 

e=^Pa (49) 

where a is the common radius of the disks. 

One of the motivations of this study was the possibility that gaseous structures attached to the 
solid wall, such as nano-bubbles, could furnish a mechanism explaining the partial slip observed by 
several investigators and it is therefore interesting to examine how the result (|49|l compares with 
available data. A full comparison would require simultaneous data for £, (3 and a. The only paper 
in which all this information is available seems to be the study by Watanabe et al. (3^, whose 
data, according to Lauga & Stone |l4j | imply a slip length of about 450 /im and an area fraction 
(3 ~ 10%. With these data, (|49J) gives a ~ 13 /xm. Rather than disks as in the present study, the 
free-slip islands in Watanabe ct al.'s work were cracks with a width of about 10 /xm and a length 
of the order 100 /im. If an equivalent radius is estimated as na 2 = 10 x 100 /jm 2 , one finds a ~ 18 
/im which is not too far from the estimate obtained from 14911 . 

The study of Simonsen et al. [24| quotes a ~ 75 nm and (3 ~ 60%. With these values, the 
estimate gives £ ~ 13 nm. Although, for such large /3's, the relation is probably not very 
accurate, this value for the slip length is in the ballpark measured by several investigators, such 
as Zhu & Granick [3(|, who report < £ < 40 nm for water, and Craig et al. [jj, who report 
< £ < 18 nm, for water-sucrose solutions. 

Wu, Zhang, Zhang, Li, Sun, Zhang, Li & Hu [34| measure a very low nano-bubblc number 
density of about 3 bubbles per 10 /im 2 , with typical radii of the order of 100 nm, which gives j3 ~ 
1% and £ ~ 0.3 nm. This is small, but not out of line with some of the existing measurements. 

The radius of surface nano-bubbles reported by Holmberg et al. is in the range 25 to 65 nm 
while that reported by Ishida et al. is of the order of 300 nm. With an area coverage of 20%, 
we can estimate a slip length between about 2 and 20 nm. Again, these numerical values are in 
the expected range. 

Tyrrell & Attard [3lJ and Steitz et al. H3] measure an area coverage of the order of 90%, which 
falls well outside the domain of applicability of our result. Unfortunately, neither group measured 
the slip length. 

Tretheway & Meinhart |28j | measured a slip length of about 1 yum, but made no estimates of area 
coverage or bubble size. With £ = 1 ^m, (|49|l gives a bubble radius a as large as 3.5 /im even for 
/3 ~ 1, and larger still for smaller (3. This is another case for which it would be of great interest to 
have some information on the surface structures. 

It is also of interest to compare our results with those of Lauga & Stone 01 obtained for flows 
in a tube with a periodic arrangement of free-slip rings perpendicular to the flow. For large tube 
radius, this arrangement should be comparable to our two-dimensional analysis. Their solution 
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is numerical, but they provide an approximate analytic expressions valid for large tube radius, 
namely 

where H is the spatial period. Upon expanding for small /?, we find 

t - y & H? (51) 

which, with the identification H/3 = a, is in precise agreement with our two-dimensional result 
(|48(l . Lauga & Stone also give a similar result for free-slip strips parallel to the flow, but this 
situation is not comparable with either one of the two that we have considered. 
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APPENDIX A: SOLUTION OF THE THREE-DIMENSIONAL PROBLEM 



We take the center of the island as the origin, with the z-axis normal to the plane and the a;-axis 
parallel to the tangential component of the traction (t) ■ n. Since the normal velocity component 
vanishes, with this choice of coordinates we require 



dv x _ s dvy 
dz dz 



0. 



(Al) 



where 



/' 



«r>-n) a 



(A2) 



Here and in the following we write v in place of (v)i for convenience. Furthermore we measure 
lengths with respect to the island radius a, although no special notation will be used to indicate 
dimcnsionlcss variables. It is convenient to adopt a system of cylindrical coordinates (r, 9, z) in 
which v x = v r cos9 — vgsin8, v y — v r sin 8+vg cos 9, in terms of which the condition (??) becomes, 
after suitable non-dimcnsionalization, 



dv r 

dz 



dv s 
dz 



S cos 6 , = — S sin 9 . 

Following Ranger [2(]| (see also 0>|2{|), we represent the velocity field in the form 



V x 



sin 9 



-x(r,z)e z + V x 



cos 9 



-ip(r,z)e, 



where e z is a unit vector normal to the plane and 



L X = LV = 



with 



dr 2 

The Cartesian velocity components follow from l|A4jl as 



ld_ 

r dr 



dz 2 ' 



v x (r,z,6) = -rd r 



1 



{dzip - X) 



cos 29 



2r 



d r (d z ip + x) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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v y {r,z,6) = -rd r 



(d z ip - x) 



sin 20 



(A8) 



v z (r, z, 6) = —d r ( -dr'tp I cos 



while, from the Stokes equation, the pressure is found as 

p(r,z,6) = fi— — Lip. 
The solution of (|A5|1 is sought in the form of Hankcl transforms with the result 



(A9) 



(A10) 



tfj = rz / e kz Ji(kr)'ij;{k)dk 
Jq 



(All) 



X = r e~ kz .h{kr)x(k)dk. 
Jo 



(A12) 



The functions ij> and \ must be determined by imposing the boundary conditions. Upon substi- 
tuting (|A11|I and (|A12|I into l|A7|l and (|A8() . we find that the no-slip condition outside the disk is 
satisfied provided that 



Ji (kr) (tp(k) + xO)) dk = - r > 1 



(A13) 



Ji(fcr) 00) - xO)) dk = r>l 



(A14) 



where d is an integration constant to be determined later. The stress condition IjAlfl inside the 
disk is satisfied provided that 



Ji(fcr) (-2ip(k) - X (k)) k dk = Sr 0<r<l 



(A15) 



Ji(fcr) (-2V>(fc) + X0))^ dk = br 0<r<l 



(A16) 



where b is another integration constant. Upon adding and subtracting, we find two pairs of dual 
integral equations for ip and x : 



Mkr)i}(k)dk = — Kr 
2r 



(A17) 



J x {kr)t})(k)dk = --{b + S)r 0< r < 1 



(A18) 



and 



r°° d 

/ J 1 (kr)x{k)dk= — Kr 
Jo 2r 



(A19) 



Ji{kr)x{k) dk = -{b - S)r 0<r<l 



(A20) 
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Both these problems have the standard Titchmarsh form 



f°° B 
I Mkr)c(k)k dk = — 1 < r 
Jo r 



(A21) 



Ji(kr)c(k) dk = Ar < r < 1 



the solution of which is (see e.g. Sneddon [2(| P-84) 



2/2 A J 5/ 2{k) B smk 



(A22) 



(A23) 



3 V 7T ^fk k 

With this result the Hankcl transforms can be evaluated in their complementary intervals finding 



J Ji(fcr)c(fc)fc dk= [b- 



AA\ 



1 



AA 



3tt J r\Jr 2 — 1 2wr 



\/ r 2 — 1 — r 2 arcsin 



r > 1 (A24) 



f 00 A A i 

/ J x (kr)~c(k)dk = rVl-r 2 +B 

Jo 3 7T 



1-VT 



< r < 1. 



The second expression is regular at r = provided that 

AA 



B 



3tt 



Upon imposing this condition on the solutions for ip and x we nn d 

d = ~S, 

so that, finally, 

$ = X 



AS h/%{k) AS k cos k - sin k 



3V2^ Vk 3tt fc 3 

The velocity field inside the disk is readily calculated from these expressions finding 

AS_ 

3?r 



v x (r,0,6) = - — Vl-r 2 v y (r t 0,9) = 0. 



The average velocity over the disk is found from direct integration: 



1 f 2lr f 1 8S* 
- / dOl rdrv x (r,Q,e) = - — 



r\J 1 — r 2 dr 



9?r 



-S 



(A25) 



(A26) 



(A27) 



(A28) 



(A29) 



(A30) 



while the y component vanishes. Although not necessary for the solution of the problem at hand, 
it may be of interest to also show explicitly the expressions for the velocity and pressure fields 
away from the disk. With the definitions: 



li 



1 



[v/(r + l) 2 + z 2 - v/(r-l) 2 + z 2 



(A31) 



h = -y{r + l) 2 + z 2 + ^(r-l) 2 + z 2 



(A32) 



the integrals can be evaluated to find (see Gradshteyn & Ryzhik ||| sections 6.621, 6.751 and 6.752) 



2Sz 2 Jii - 1 „„ 45 f I ~ . ( 1 

v x (r, z, 9) = -_r ( ^_^ cos2fl - - 1 - 1\ - zarcsm ^_ 



2Sz y^T~T . /I 

-TTn tt ^ arcsm — 

3tt l(£i-€ 2 ) 



(A33) 
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v v (r,z,6) = -^ Jjjgl^ rinM (A34) 



„{r,z,6) = (- if^f 1 ) cos20 (A35) 



3tt \ (i 2 2 -£l)r 



^ = ^[-W^w)™ 26 (A36) 



APPENDIX B: SOLUTION OF THE TWO-DIMENSIONAL PROBLEM 

In this case it is convenient to adopt as fundamental length and a Cartesian system of 
coordinates with x along the plane direction and y along the normal. The Boundary conditions of 
the Stokes problem for v x and v y become 

v x (x,0) = \x\>l (Bl) 

d y v x (x,0) = S \x\<l. (B2) 

v y (x, 0) = — oo < £ < oo . (B3) 

where 

S = -~(r X y). (B4) 

We introduce a stream function ip in terms of which 

v x (x,y) = dyip, v y (x,y) = -d x ip (B5) 

and 

oj = d y v x - d x v v = Aip. (B6) 
The vorticity u> is harmonic and can be written as a Fourier integral in the form 

poo 

w(x,y) = dkexp(ikx)uj{k) e -^ y . (B7) 



By introducing the Fourier transform if)(k,y) of the stream function, substituting into (|B6|I . and 
integrating, we hnd 

fe) = -^e-l*l» (B8) 



after elimination of an integration constant on the basis of (IB3|) . With this result, the boundary 
condition (|B1|1 becomes 

/oo 
dkexp(ikx)u(k) = S \x\ < 1 (B9) 
-CO 



and 1B2 



dk exp{ikx) = \x\ > 1 . (BIO) 
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Upon writing (|B9() for x and -x and adding or subtracting, we find 

/oo 
dk cos(kx)Cu(k) = 5 < x < 1 (Bll) 
-OO 

/OO 
dfcsin(fcz)£(fc) = 0<x<l. (B12) 
-oo 



Proceeding in a similar way with (|B10|) we have 

dkcoa(kx)^r = Kx (B13) 
l«l 

dksm(kx)^$- = 1< x. (B14) 
\ k \ 

If in (|B11|) we separate the integration range into — oo < k < and < k < oo we find 

d/s cos(kx)Cu + = S < x < 1 <!>+ = cD(fc) + fc) (B15) 
whereas (|B12|) gives 

dfc sin(fcz;)(I>_ = < x < 1 cD_ = w(fc) — ui(— k). (B16) 







Similarly 



(i 



dkcos(kx)— = 1 < x (B17) 
fc 



f 00 Co 

/ dfcsin(fcc) — = Kx. (B18) 
Jo fc 

Since the problem for £>_ is completely homogeneous, this quantity must vanish so that Cb{k) is 
even in k and, therefore, real. We are thus led to the pair of dual integral equations 



,00 1 

/ dk cos(kx)uj = —S < x < 1 
7o 2 



(B19) 



dkcos{kx)y =0 1 < x. (B20) 



This is a standard problem with the solution (see e.g. Sneddon j2g p. 84) 



Co{k) = -SJi(k) (B21) 

from which the velocity on the boundary follows as 

g g 

v x (x, 0) = — — cos (arcsinx) = — — yl — x 2 x < 1 (B22) 

so that 

±J v x (x,0)dx = ~S. (B23) 
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This result coincides with that derived by different means in Philip [jj]]- As before, it may be of 
some interest to show the explicit results for the velocity and pressure fields. One has 



Vx(x,y) 



~j dk cos{kx) fjij -y^j u{k)e~^ y (B24) 



S f °° 

v y (x,y) = -- I dk sin(kx)yCj(k)e^ klv (B25) 



p(x,y) = -25 / dk sin{kx)Co{k) e- kv . (B26) 
Jo 

The integrals can be evaluated to find 

v x (x,y) = -^R( x ,y)-^d y R(x,y) (B27) 



v v (x,y) = ^-d y I(x,y) (B28) 



p(x,y) = Sd y R(x,y) (B29) 



with 



2 



/ (1 - x 2 + y 2 ) + ^(1 - x 2 + y 2 ) 2 + 4x 2 j/' 

R{x,y) = -y+\ — ~ ( B3 °) 



/ -(l-.T 2 +y 2 ) + ^(l-^+^^+^V 

7(a;,2/)=x+l/ ^— . (B31) 
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